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ABSTRACT 

The rich structure of the A 1 !! — X 1 'E + system of carbon monoxide accounts for many of the spectral 
features seen in long slit i?ST-STIS observations of comets 153P/Ikeya-Zhang, C/2001 Q4 (NEAT), 
and C/2000 WMj (LINEAR), as well as in the H ST-GHRS spectrum of comet C/1996 B2 (Hyaku- 
take). A detailed CO fluorescence model is developed to derive the CO abundances in these comets 
by simultaneously fitting all of the observed A — X bands. The model includes the latest values for the 
oscillator strengths and state parameters, and accounts for optical depth effects due to line overlap 
and self-absorption. A complete fitting of the CO spectral features using this model leads to the first 
identification of a molecular hydrogen line pumped by solar H I Lyman-/? longward of 1200 A in the 
spectrum of comet 153P/Ikeya-Zhang. Pumping by strong solar lines also plays an important role in 
CO fluorescence, as shown by the detection in the spectrum of comet C/1996 B2 (Hyakutakc) of bands 
from the (14-v") and (9-v") progressions pumped by solar H I Lyman-o and O I A1302, respectively. 
Using spectra extracted at increasing distances from the comet nucleus or averaging over increas- 
ing effective apertures, the model fits yield radial profiles of CO column density that are consistent 
with a predominantly native source for all the comets observed by STIS. The derived CO production 
rates are 1.54 ± 0.09 x 10 28 molecules s _1 for 153P/Ikeya-Zhang, 1.76 ± 0.16 x 10 28 molecules s _1 
for C/2001 Q4 (NEAT), and 3.56 ± 0.20 x 10 26 molecules s" 1 for C/2000 WMi (LINEAR). In the 
absence of spatial information for comet C/1996 B2 (Hyakutake), we estimate a CO production rate 
of 4.97 ± 0.07 x 10 28 molecules s"\ assuming an entirely native source. The CO abundances relative 
to water in these comets span a wide range, from 0.44 ±0.03% for C/2000 WMi (LINEAR), 7.2 ±0.4% 
for 153P/Ikeya-Zhang, 8.8 ±0.8% for C/2001 Q4 (NEAT) to 20.9 ±0.3% for C/1996 B2 (Hyakutake). 
For comets C/2000 WMi and C/2001 Q4 we can compare these results with those derived from nearly 
simultaneous observations by the Far Ultraviolet Spectroscopic Explorer. 

Subject headings: cometsundividual (153P/Ikeya-Zhang, C/2001 Q4 (NEAT), C/1996 B2 (Hyakutake), 
C/2000 WMi (LINEAR)) — carbon monoxide — ultraviolet: solar system 



1. INTRODUCTION 

Cometary ices have undergone little or no processing 
since their formation in the solar nebula and thus they 
represent an important clue in understanding the condi- 
tions in the early Solar System. However, this interpre- 
tation is complicated by the radial mixings in the pro- 
toplanetary disk and the dynamical evolution of comets, 
expelled to the outer Solar System from their original 
formation site. No pattern for the cometary composi- 
tion and activity has eme r ged from the o bjects studied 
to date (jBiver et al.ll2002t ICrovisied[2007| ) and a larger 
sample is needed. The composition of cometary ices has 
been studied through spectroscopic observations or in 
situ measurements. In the analysis of cometary spectra 
it is important to establish whether the observed chemi- 
cal species originate in the nucleus (native source) or in 
the coma (extended source). For this purpose, spatially- 
Electronic address: roxana@pha.jhu.edu 



resolved observations are needed to derive the radial com- 
position of the gas outflow. 

Carbon monoxide (CO) is one of the most abun- 
dant molecules, whose mixing rati o relative to wa- 
ter varies greatly among comets ([Biver eTaTl 12001 
iBockelee-Morvan et al.ll2004l ). As a highly volatile com- 
pound, with a sublimation temperature of 24 K, CO ice is 
thought to be a sensitive tracer of the temperature in the 
environment in which the comets formed. However, its 
origin as a native compound or a daughter product, and 
the correlation between the two sources is still not un- 
derstood. Infrared observations of the CO ro- vibrational 
lines conclude that the CO sou rce was mainly native in 
comet C/1996 B2 (Hyakutake) (jDiSanti et al.ll2003l ) and 
significantly extende d for comet C/1995 01 (Hale-Bopp) 
(jDisanti et alJ 12001 s ). This paper complements previous 
findings with the first CO radial profiles constructed from 
ultraviolet spectra. Care must be taken in interpreting 
the brightness profiles, since saturation close to the nu- 



2 



Lupu et al. 



cleus could mimic the presence of an extended source. A 
comprehensive fluorescence model is developed to derive 
reliable CO column densities and estimate CO produc- 
tion rates from cometary spectra obtained by the Space 
Telescope Imaging Spectrograph (STIS) and the God- 
dard High Resolution Spectrograph (GHRS) on the Hub- 
ble Space Telescope (HST). Using the long-slit capabili- 
ties of the STIS instrument, this paper offers evidence 
that the native source of CO is dominant in three long 
period comets. 

CO ultraviolet fluorescence in cometary spectra was 
first detected during sou nding rocket observatio ns of 
comet West (C/1975 VI, iFeldman fc B7u^fl97fih and 



has b een subs equently obser ved by IU E (ITozzi et al 



1998), HST (jWeaverl [l998h . FUSE (|Feldman et al 



20021) . the Hopkins Ult raviolet Telescope o n the Astro- 1 



Space Shuttle mission (IFeldman et al.lll99lD, and rockets 
(IWoods et all 1198ft ISahnow et all 119931: iMcPhate et all 
I1999T ) . The most important spectral features of CO in 
the ultraviolet (UV) are its electronic transitions belong- 
ing to the A — X, B — X and C — X systems, or to 
the forbidden Cameron bands. Fluorescence is the main 
emission mechanism for the A 1 !^ — X 1 Y 1 + system (1300 - 
1900 A), C 1 E+ - X X Y,+ system (0-0 band at 1087.9 A), 
and - system (0-0 band at 1150.5 A). The 

Cameron bands, a 3 H.-X 1 E + (1900- 2800 A), are mainly 
excited by electron impact and photodissociation of CO2 
(jWeaver et aljl 994). Although observations of the A— X 
or Fourth Positive Group of CO in the UV spectra of 
comets have a long history, their interpretation has been 
difficult compared to that of the B—X and C— X systems 
at shorter wavelengths that have been observed more re- 
cently. 

The spatial information offered by the high resolu- 
tion long slit 7JST-STIS spectra of comets 153P/Ikeya- 
Zhang (C/2002 CI), C/2001 Q4 (NEAT), and 
C/2000 WMi (LINEAR) shows that close to the nucleus 
the CO emission in the A — X system is self-absorbed. 
This follows from the observed change in the relative in- 
tensities in various vibrational progressions as the offset 
from the comet center increases (see § 3). Self-absorption 
makes it difficult to derive a reliable value for the CO 
column density in the absence of a model that takes into 
account optical depth effects. 

A detailed model of the A — X system is needed 
to track the effects of saturation and self-absorption. 
We constructed a database containing ~ 10 5 transi- 
tions between the first 50 rotational levels of each of 
the 37 vibrational levels of X 1 T, + and the 23 vibra- 
tional levels of A 1 ]!, taking into account the energy 
shifts and mixings of the transition probabilities due 
to interactions between the different parity sublevels 
(|Le Floch et all Il98ft iMorton fe Noreaul I1989D . Using 
a sim ple approximation for fluorescence in subordinate 
lines ([Liu fc Dalgarnol ll996f). expanded with a compre- 
hensive treatment of self-absorption, our model offers 
an excellent fit to the data. The A — X fluorescence 
model is used to derive spatial profiles of the CO col- 
umn density for the three comets observed by STIS. 
We find that the column density profiles are consis- 
tent with a dominant native source in all comets. The 
resulting production rates range between 3.56 x 10 26 
and 1.76 x 10 28 molecules s^ 1 and are corroborated 



by the results from high resolution FUSE observations 
of C/2001 Q4 (NE AT) and C/2000 WMj (LINEAR) 
(I Weaver et al.ll200l . The HST-GBRS observations of 
comet C/1996 B2 (Hya kutake), a comet with a high CO 
production rate (jBiver et all 119991 : iDiSanti et aLll2003l ) 
and thus strongly affected by large optical depth effects, 
require the addition of geometrical corrections to the sim- 
ple plane-parallel atmosphere model in order to repro- 
duce the relative line strengths. 

The HST observations are described in § 2. Details 
about the model and its simplifying assumptions are 
found in § 3. Detailed data analysis follows in § 4, and 
a discussion of the results is given in § 5. We conclude 
with a summary in § 6. 

2. OBSERVATIONS 

The HST observations are summarized in Table [TJ 
Comets Ikeya-Zhang and C/2000 WMi (LINEAR) 
were observed near times of high solar activity, while 
C/2001 Q4 (NEAT) and Hyakutake were observed closer 
to solar minimum. The -ffST-STIS observations used 
the G140L grating and the 52X0.2 aperture (25" x 0'.'2 
for L-mode MAMA), resulting in a spatial resolution of 
~0'/l and a spectral resolution of ~4 A in the wavelength 
range 1150-17 30 A. The STIS ins trum ent performance 
is desc ribed in lKimble et al | (1 19981) andlWoodgate et aTl 
(fl998l) . The HST-GBRS (|Heariet^IJ1995|) observations 
of comet C/1996 B2 (Hyakutake), covering the 1290- 
1590 A bandpass with a spectral resolution of ^4 A, 
do not provide spatial information. The GHRS Large 
Science Aperture (LSA) translates into a l'/74 x 1'.'74 
post-COSTAR projected area. Given that FUSE and 
HST-STIS observations of comets C/2000 WMi (LIN- 
EAR) and C/2001 Q4 (NEAT) were nearly simulta- 
neous, we choose to include in Table D] a summary 
of FUSE observations fo r completeness ()Feldman et all 
l200llWeaVer et alJl2002l ). 



3. FLUORESCENCE MODEL FOR THE A 1 n 
SYSTEM OF CO 



The Fourth Positive Group of CO, A 1 !! - A 1 E+, 
has non-negligible overlap integrals for most of the non- 
diagonal vibrational transitions. While for the C — X 
and B — X systems it is enough to model just the (0-0) 
and (1-0) bands, for the A — X system we must take 
into account all bands connecting all 37 vibrational lev- 
els of the X 1 T, + state with all 23 vibrational levels of 
the A 1 ]! state. Using only the first 50 rotational lev- 
els, which should be sufficient for typical physical con- 
ditions in cometary comae, the final database contains 
almost 10 5 transitions . The l atest values for the pa- 
rame t ers of the A 1 !! (jKuruczl 119761 : IMorton fc Noreaul 
1989 1 : iLe Flochl fl99lT ) and X X S+ states (|George et all 



19941 ) were used to derive th e energy levels and transi- 



tion wavelengths. Following IMorton fc Noreaul (|1989l ) 
the rotational transition probabilities Ay u were ob- 
tained from the band transition probabilities A 1 „ 
(iBeegle et all 119991: Hidelsberg et all 119991 : iBorges et ail 
12001b iKirbv fc Cooper! Il989f ). Even though there are 
three rotational branches (P, Q and R, AJ = 0, ±1) pos- 
sible from the same rotational level of the A 1 !! state (J'), 
due to the parity selection rule, the excitation rates and 
branching ratios for the P and R branches do not mix 
with those for the Q branch. Each J' level of the ^4 1 n 
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TABLE 1 
Comet observations. 



Dataset 


Date & Time 


Exposure Time 






A 


Instrument 


Mode 




(UT) 


(a) 


(AU) 


(km s _1 ) 


(AU) 






C/1996 B2 (Hyakutake) 














Z35FN602T 


1996-04-01 13:09:00 


1305.600 


0.885 


-38.51 


0.259 


iTST-GHRS 


G140L;2.0 


Z35FN604T 


1996-04-01 13:39:00 


zl f .oUU 


U.ooo 


— oo.ol 


U.ZOU 


no 1 -Kjtiixo 


(~* 1 A HT -O n 

1_j14U-L,Z.U 


Z35FN702T 


1996-04-01 14:44:00 


lOUJ.UUU 


U.OO'i 


oo.oo 


262 


TJ C?rp pTTDC 
11 O 1 -\J LLIXO 




Z35FN704T 


1996-04-01 15:14:00 


917 finn 

ZL ( .uuu 


u.ooo 


oo.oo 


262 


TJ qrp pTTDC 

11 O 1 -KjrLIXO 


n 1 aot -9 n 


C/2000 WMi (LINEAR) 














B0500301000 


2001-12-07 08:50:00 














to B0502401000 


2001-12-10 06:48:00 


36,467 


1.120 


-28.30 


0.340 


FUSE 


LWRS 


O6GR12010 


2001-12-09 21:33:00 


1800.194 


1.085 


-28.26 


0.357 


HST-STIS 


G140L;52x0.2 


O6GR03010 


2001-12-09 23:09:10 


1440.197 


1.084 


-28.26 


0.358 


HST-STIS 


G140L;52x0.2 


O6GR11010 


2001-12-10 00:45:00 


1800.198 


1.083 


-28.26 


0.359 


HST-STIS 


G140L;52x0.2 


153P/Ikeya-Zhanj 
















O8FY01010 


2002-04-20 07:28:00 


1800.199 


0.887 


29.07 


0.426 


HST-STIS 


G140L;52x0.2 


O8FY02010 


2002-04-20 10:40:00 


1800.197 


0.889 


29.08 


0.426 


HST-STIS 


G140L;52x0.2 


O8FY03010 


2002-04-21 07:30:00 


1440.197 


0.904 


29.15 


0.422 


HST-STIS 


G140L;52x0.2 



C/2001 Q4 (NEAT) 

E1390101000 2004-04-24 00:39:00 

to E1390501000 2004-04-24 23:09:00 68,282 1.030 -10.80 0.510 FUSE LWRS 

O8VK04010 2004-04-25 20:03:00 1800.199 1.024 -10.26 0.473 HST-STIS G140L;52x0.2 

O8VK01010 2004-04-26 00:49:00 1683.008 1.023 -10.18 0.468 HST-STIS G140L;52x0.2 

O8VK07010 2004-04-30 00:51:00 1800.200 1.002 -8.37 0.386 iTST-STIS G140L;52x0.2 



state is split into 2 opposite parity sublevels (A-doubling) 
that have different energies and inter actions with neigh- 
borin g levels of other electronic states (|Morton fc Noreaul 
1989). The energy level shifts and the changes in the life- 
times due to the se interactions were taken into account 
when available dMorton fc Noreaul Il989t iKittrell et al.l 
[19981 : iLe Floch et al.lll987f T 

In a first approximation, the fluorescent emission in the 
Fourth Positive G roup of CO is modeled following the 
prescription from iLhi fc Dalgarnol ([1996). This model 
successfully accounts for saturation in the absorption of 
the exciting radiation, using a Voigt profile with line 
overlap. A quick look at the data, however, shows that 
this approach only partially accounts for optical depth 
effects, and we must extend the model to include self- 
absorption in the fluorescent cascade. The STIS spectra 
of comet 153P/Ikeya- Zhang shown in Figure Q] are aver- 
aged over effective apertures of increasing size, centered 
on the comet nucleus. Self-absorption is easily recog- 
nized by noting that the bands connecting to the ground 
vibrational level (0-0, 1-0, 2-0, 3-0, 4-0, 5-0) show lit- 
tle decrease in brightness as we increase the integrated 
slit area and the average column density becomes lower, 
while the other bands from the progressions originating 
in the same upper level (see the unblended 0-1, 1-3, 2- 
2, 2-3, 4-1 and 5-1) decrease strongly in intensity, so 
that at low column densities the relative line ratios are 
in agreement with the optically thin limit. 

The treatment of self-absorption is made under the as- 
sumption of local thermodynamic equilibrium (LTE) , us- 
ing the fact that the photon mean free path corresponds 
to an optical depth t of one. Given that the gas tem- 
perature in the coma is less than 100 K for our obser- 
vations, under LTE conditions we need only to take into 
account absorption out of the lowest (0) vibrational level 
of X 1 Yi + . Unlike for H2 for example, the ro- vibrational 
levels of the ground state of CO have short lifetimes, 
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Fig. 1. — Spectra of comet 153P/Ikeya-Zhang extracted from 
apertures that are 1, 2, 4 and 6 arcseconds wide in the spatial 
direction, top to bottom, centered on the comet nucleus, showing 
the effects of self-absorption (see text). The strong features near 
1564 A and 1657 A are emissions from atomic carbon. 



making LTE a good approximation. The lines connect- 
ing to v" = for which r at line center is greater than 
one are considered self-absorbed. Only an effective col- 
umn density Nq = l/c^„ e _ center , smaller than the to- 
tal absorbing column N, will contribute to the observed 
emission in a self-absorbed transition i — j. The excess 
brightness, due to absorption by the total column den- 
sity N, is then redistributed among the optically thin 
lines originating in the same upper level i according to 
the branching ratios, as seen in Figure [5J 

The model is constructed in the approximation of a 
uniform plane-parallel atmosphere. This approximation 
is of limited use, due to the spherical symmetry of the 
comet atmosphere and the Sun-comet-Earth geometry. 
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Fig. 2. — Illustration of the CO fluorescence model. The absorp- 
tion profile for the A— X (1-0) band is applied to the solar spectrum 
(upper panel) and the subsequent emission in the (1 — v") progres- 
sion is shown before and after including self-absorption (middle 
and lower panels, respectively). 
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The breakdown is more apparent for larger column densi- 
ties, such as in the case of comet C/1996 B2 (Hyakutake), 
discussed in § 4.2. In the optically thin regime, valid at 
larger distances from the comet nucleus, as well as for un- 
resolved objects, such geometrical effects are negligible. 
Introducing geometrical corrections in our model can be 
done by allowing the column density entering the absorp- 
tion step to be different from the column density used for 
emission and self-absorption. This procedure accounts 
for the fact that the projected column density along the 
line-of-sight is larger than the column density towards 
the exciting radiation source. This difference leads both 
to the decrease in line saturation and to the enhancement 
of the optically thin lines versus the optically thick ones 
due to more self- absorption. The method is very robust, 
providing line intensities in agreement with the data. 

4. DATA ANALYSIS 

Although the brightness of the A — X bands for the 
same comet should differ slightly from one observation to 
another due to varying comet heliocentric and geocentric 
distances, as well as due to possible periodic variations in 
the volatile vaporization rate, this variation is within the 
error bars of the observations and the use of an averaged 
STIS spectrum for each comet in order to improve the 
signal-to-noise ratio is warranted. We also select only 
the datasets that do not show significant deviations in 
background and intensity from one another. 

High resolution solar spectra from the Ultraviolet Spec- 
trometer Pola rimeter Experiment on the S olar Maxi- 
mum Mission (|Tandberg-Hanssen et al.|[l981h were used, 
scaled to match the solar activity at the time of comet 
observations as deduced from UARS /Solstice solar flux 
measurements (]Rottman et al.l l200I). The solar spec- 
trum is shifted according to the comet motio n relative 
to the Sun (Swings effect. iDvmond et al.l H989). For the 
so lar H I Lyman-q and -0 lines the SOHO/SVMER data 
of iLemaire et al.l (|2002f ) were used. 

The only free parameter of the model is the column 
density, which is varied over a grid of values until the 
best fit is found. All model parameters used for each 
comet are summarized in Table [2j The values for the 
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Fig. 3.— HST-ST1S spectral image of comet 153P/Ikeya-Zhang, 
showing the line brightness variation along the slit. The spectral 
region containing the strong gcocoronal H I A1216 and O I A1302 
has been excluded. The carbon multiplets at 1561.0 A and 1657.6 A 
arc relatively constant due to the extended source component. The 
zero point in the vertical direction marks the location of the bright- 
ness peak, chosen as the center of our integration bins. 



rotational temperature and Doppler b parameter were 
chosen in agreement with infrared and radio measure- 
ments, when available. The Doppler b parameter is given 



by b = ,Jvf herm + v 2 outflow a Voutfiow, where v outflow 

is the source of non-thermal line broadening, and the 
thermal velocities are comparable to the uncertainties in 
v out flow Molecular lines of water and other molecules 
are well resolved by radio observations (usually to better 
than 0.1 km s _1 ), and have a Ga ussian profile, reflecting 
the symmetric outflow of the gas dLecacheux et ai1l2003t 
IWouterloot et alJll998t iBiver et al.lll999h . Our"observa- 
tions have too low a resolution to be used to determine 
the Doppler line widths directly, so we rely on the pub- 
lished values for the expansion velocity of the gas, de- 
rived from the line profiles after correcting for thermal 
and instrumental broadening. 

4.1. HST-STIS Observations 

The high spatial resolution of the iTST-STIS instru- 
ment is illustrated by the spectrum of comet I53P/Ikeya- 
Zhang in Figured The CO A — X lines form the major- 
ity of the observed features. Their intensity decreases 
rapidly along the slit, in contrast with the extended 
C I A156I.0 and AI657.6 multiplets. Using the spatial 
information available, we derive the CO spatial profile 
for each comet by fitting spectra extracted from regions 
of l'/5 width at increasing offsets from the comet cen- 
ter. The selected regions sample areas of varying column 
density, for which our model estimates an average value. 
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TABLE 2 

Model parameters. 



Comet Name 


Trot 


Doppler b parameter 


Solar Activity 


N C o Range 




(K) 


(km s^ 1 ) 




(10 14 cm" 2 ) 


153P/Ikeya-Zhang 


82 a 


0.91 b 


max 


61.3-1.49 


C/2001 Q4 (NEAT) 


68 c 


0.79 d 


min 


68.4-1.86 


C/2000 WMi (LINEAR) 


77 c 


0.72 b 


max 


0.935-0.327 


C/1996 B2 (Hyakutake) 


72° 


2.0 f 


min 


145 



a From 74xr-°- 93 K dependence. IDello Russo et al.l [120041) . b IBiver et al.l j2006h . c Temper- 
ature of cold component, FUSE observations (see text). d Estimated from outflow velocity 
0.8xrv a5 . c From 63xr~ 1M K dependence, IDiSanti et"aLl [|2003T ). f Within the range of 



values given bvlWoutcrloot et al. (1998). 

The comet nucleus is identified with the center of bright- 
ness, located at zero in Figure [3] The innermost region 
extracted, centered on the nucleus, is noisier due to the 
small integrated area. For intermediate regions the signal 
is better, due to the averaging of two regions, symmetric 
about the comet center. For each region, the background 
subtraction is performed by fitting a quadratic polyno- 
mial to selected points from feature-free intervals. These 
points were selected such that the outliers are discarded 
and the resulting polynomial fit is optimal. 

The best-fit column density for each region and its 
standard deviation are derived by minimizing the x 
statistics, taking into account the errors in background 
subtraction. The range of column densities obtained for 
each comet observed by STIS is listed in Table[2j together 
with the average value for comet C/1996 B2 (Hyaku- 
take). The results are further compared with the Haser 
native source model ()Haserl 119571 lOpal fc Carruthersl 
119771) . with an outflow velocity of 0.85 r~ 5 km s" 1 , 
where is the comet-S un distance in AU ( Budzic n~et al.l 
[1991 IBiver elallfl999l) . and a CO lifetime of 2 x 10 5 s. 
We derive the CO production rate for each comet by a 
least squares fit of the Haser model to the radial column 
density profile. The native source model is integrated 
over rectangular regions matching the l'/5 spectral ex- 
tractions along the STIS slit. The resulting production 
rates and their magnitude relative to water are listed in 
Tabled 

4.1.1. 153P/Ikeya-Zhang 

The values for the rotational temperature and 
Doppler b parameter are 82 K and 0.91 km s _1 
respectively , derived from infrare d and radio mea - 
surements (jDello Russo et al.l 12004 IBiver et al.l I2006D . 
Spectra extracted from l'/5 intervals were fitted us- 
ing CO column densities ranging from 6.13xl0 15 to 
1.49 xlO 14 cm -2 , as shown in Figure |U The data shown 
are obtained by averaging the first and third STIS obser- 
vations (Table [J) , and the best fit model is overplotted 
in red. The second observation was not included due to 
the background mismatch with the other two. The de- 
rived values for the CO column density in each region are 
represented by stars in Figure [5j plotted as a function of 
the distance from the comet nucleus. The error bars are 
given by the la confidence level from the \ 2 statistics. 
Fitting to this radial profile a native source model inte- 
grated over rectangula r regions with the same cov erage as 
the extracted spectra (|Qpal fc Carrutherslll977i ). we ob- 
tain a production rate of 1.54± 0.09 x 10 28 molecules s — 1 . 



The native source model for the derived production rate 
is shown as a continuous line in the same figure. The re- 
sulting CO production rate relative to water is about 
7.2 ± 0.4%. The water production rate (2.15 x 10 29 
molecules s _1 ) was obtained from a vectorial model fit 
to an iJST-STIS observation of the OH (0-0) band made 
on 2002 April 21 at 12:19 UT. 

The residuals after subtracting the fluorescence spec- 
trum for the CO A — X system reflect the contributions 
of atomic species and allow the first detection of H2 at 
wavelengths longward of 1200 A. The H2 spectrum con- 
sists of the P(l) lines of the Lyman (6-v") progression 
pumped by solar H I Lyman-/?. The H2 lines with v"= 
1-3 were first detected in the FUSE spectru m of comet 
C/2001 A2 (LINEAR) (jFeldman et alJ 12001 ). The lines 
at longer wavelengths, including the strongest one in the 
progression (v"= 13 at 1607.5 A), remained undetected 
due to the abundance of CO features and lower resolu- 
tion of the STIS instrument. The solar maximum value 
of flux, together with a velocity shift that placed the (6- 
0) P(l) line at the peak of the line, made it particularly 
fortuitous to detect the (6-13) P(l) line in comet Ikeya- 
Zhang. Figure [5] shows the shape and intensity of the 
solar Ly/3 at minimum and maximum activity and the 
Doppler shifts of the H2 absorption line corresponding 
to each comet. Although H 2 l ines pumped by Ly/3 were 
detected in comet C/2001 Q4 (jFeldman et al.ll2004D . due 
to the large negative Doppler shift and the low solar ac- 
tivity, in the STIS bandpass the signal-to-noise for the 
H2 lines is too low to warrant a detection. A comparison 
of the residuals for the two comets after subtracting the 
CO fluorescence model is shown in Figure [7] The STIS 
spectrum was integrated over a 4" wide region centered 
on the comet nucleus, and the CO column density for the 
two comets (N = 3.07 x 10 15 cm" 2 and 3.49 x 10 15 cm" 2 , 
respectively) has been obtained through the x 2 mini- 
mization procedure. The residuals for comets C/2000 
WMI and Hyakutake (not included in the figure) show 
no evidence for H2, which can be explained by the large 
negative Doppler shifts with respect to the solar Ly/3, the 
low outgassing rate of C/2000 WMI and the minimum 
solar activity at the time of Hyakutake observations. 

The synthetic H2 fluorescence spectra shown in red 
in Figure [7] are constructed under the assumption of 
an H2O photodissociation source for H2. According 
to the dissociation model, H2 is rotationally hot (100- 
300 K) and its rotational levels are in statistical equi- 
librium, with an ortho/para ratio of 3, similar to water 
(jBudzien et al.lll994t Ivan Harrevelt k van Hemertj|2000l: 
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Fig. 4. — HST-STIS spectra of comet 153P/Ikeya-Zhang extracted from 1.5" intervals at increasing offsets from the center of brightness, 
assumed to be the location of the nucleus. The red line represents the model spectrum of the CO A — X system for the best fit column 
density for each offset. The corresponding offsets and the best fit column densities are indicated on each panel. Other emission features 
belong to C II, C I, O I, S I, and H2, as indicated in Figure 171 



iBonev et al.ll2007t ). The uncertainty in the exact value 
for the rotational temperature has little impact on the 
emerging spectrum, as the population of the J=l level 
of the ground state does not differ significantly for ro- 
tational temperatures ranging from 100 to 300 K. The 
fluorescence efficiencies, or g-fact ors, used in the mode l 
were revised from those cited bv iFeldman et al.l (2002), 
using the sola r Lyman-/3 profiles ob tained over the past 
solar cycle bv iLemaire et al.l ()2002[ ). and accounting for 
the comet's heliocentric velocity. The quality of the data 
does not warrant a x 2 fit f° r t ne H2 fluorescence model, 



so we restrict ourselves to making rough estimates. Using 
an H2 column density of 1.0 x 10 14 cm -2 at a rotational 
temperature of 200 K we obtain a reasonable agreement 
with the residuals for 153P/Ikeya-Zhang (Figure[71 upper 
panel). As was the case for comet C/2001 A2, this value 
for the H2 column is consistent, within the rather large 
uncertainties in both the data and the models, with an 
H2O photodissociation source. Constraining the H2 pro- 
duction helps our understanding of this H2O dissciation 
channel, for which little laboratory data is available. 
Under our choice of parameters, the model for 
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Fig. 5. — Spatial profile of the CO column density for comet 
153P/Ikeya-Zhang. The stars with error bars represent the column 
density values derived from averaged spectra using the method 
presented in the text. The continuous line represents the column 
density profile predicted by the native source model when aver- 
aged over slits covering the same areas as the observations, with a 
production rate of 1.54xl0 28 molecules s _1 . 
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Fig. 6. — Profile of solar Ly/3 at maximum (August 2001) and 
minimum solar activity (July 1996). The velocity shifts of the 
absorbing H2 (6-0) PI line are shown for each comet. 



153P/Ikeya-Zhang predicts about 5.5 R for the v"= 7,9 
and 11 lines, at the level of the errors due to noise and CO 
subtraction, leaving only the (6-13) P(l) line detectable 
at a ~3cr level, with 18.4 R. This line is not detected in 
any other comet from our sample. For comparison, we 
model the H2 fluorescence for comet C/2001 Q4. We use 
the same column density of H2, as both C/2001 Q4 and 
153P/Ikeya-Zhang have similar water production rates 
(see notes to Table [3]). Due to the low solar activity and 
the unfavorable Doppler shift, the predicted line intensi- 
ties are too low to be detected, consistent with the resid- 
uals (Figure [71 lower panel). The discrepancies between 
the H2 models and the residuals are due to the uncertain- 
ties in the CO fluorescence model itself as well as to the 
noise in the data increasing towards longer wavelengths. 

The strongest atomic lines seen in the residuals are 
also indicated in Figure [Jj In addition to t hose lines 
usually seen in comets (jMcPhate et al.lll999h . we also 
note the presence of the S I ( 1 D — 1 D°) transition at 
1666.7 A. This transition is analogous to O I ( 1 D — 1 D°) 
at 1152.2 A (|Feldman et al ll2002l ) and C I ( X D - 1 P°) at 
1930.9 A (jTozzi et al.lll998h . Using a g-factor of 1.92 x 
10 -5 photons atom" 1 s" 1 at 1 AU, the ~30 rayleigh 
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Fig. 7. — Residuals from the spectrum of 153P/Ikeya-Zhang and 
C/2001 Q4 after subtracting the best fit CO model (JV = 3.07 X 
10 15 cm -2 and 3.49 X 10 15 cm -2 , respectively). The comet spectra 
were extracted from a region of 4" width, centered on the comet 
nucleus. The red line is the predicted H2 fluorescence spectrum 
pumped by solar Lyman-/3 for a column density of 1.0 X 10 cm -2 
and a rotational temperature of 200 K. Other atomic contributions 
are shown in blue. 

brightness corresponds to an average column density of 
SI 1 !? atoms in the aperture of 1.24 x 10 12 cm" 2 . Since 
the lifetime of the metastable 1 D state of sulfur is only 
28 s, this requires that within the aperture l D atoms be 
produced at a rate of 1.9 x 10 26 s" 1 (0.09% relative to 
water). Collisional de-excitation of X D near the nucleus 
would raise this number. The likely source of these atoms 
is the photodissociation of sulfur-bearing molecules such 
as H2S or CS2, which must be produced at a rate greater 
than 0.1% relative to water. 

4.1.2. C/2001 Q4 (NEAT) 

Adjusting the model parameters to the conditions of 
NEAT Q4 observations we derive in the manner de- 
scribed for 153P/Ikeya-Zhang a CO production rate of 
1.76 ± 0.16 x 10 28 molecules s" 1 , or 8.8 ± 0.8% relative 
to water. For the water production rate, we derive a 
value of 2.0 x 10 29 molecules s" 1 from STIS observations 
on 2004-04-23 21:39 UT. The CO model was fitted to 
the average of the first two STIS observations (Table [T]). 
For the third observation the detector background level 
is much higher than for the other two. The adopted 
values for the rotational temperature and Doppler b pa- 
rameter are listed in Table [2] FU SE observations of 
the CO C - X Hopfield-Birge system at 1088 A reveal 
a band profil e consistent with a tw o component temper- 
ature model dFeldman et al.ll2002f h The hot component 
(^600 K) is believed to describe an extended CO sourc e 
due to the dissociation of CO2 (Feldma n et al.l 12006b). 
We use for our model the temperature of the cold com- 
ponent, estimated at ^68 K, characteristic of the native 
CO source which dominates at the smaller cometocentric 
distances probed by STIS. Lacking direct radio measure- 
ments of the line widths, we use a value of 0.79 km s" 1 for 
the Doppler b parameter, based on the outflow velocity 
0.8 xr^ . The radial profile of the best-fit column den- 
sities and the native source model are plotted in Figure[5] 
as stars and continuous line, respectively. 
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Fig. 8.— Same as Figure [5] for comet C/2001 Q4 (NEAT). 
The production rate for the native source model was 1.76 X 
10 28 molecules s — 1 . 

Using a two-com ponent fit to the CO C — X band ob- 
served by FUSE (jFeldman et all I2002D we derive CO 
column densities of 1.2 ± 0.3 x 10 14 cm" 2 and 2.38 ± 
0.08 x 10 13 cm~ 2 for the cold and hot component, re- 
spectively, averaged over the entire 30" x 30"slit. A na- 
tive source for the cold component requires a production 
rate of 1.36 ± 0.40 x 10 28 molecules s" 1 , or 6.8 ± 2.0% 
relative to water. The solar flux pumping the C — X 
fluorescence was base d on quiet sun w hole disk data 
from SOHO/SVMER (|Curdt et al.ll2001l ). normalized, at 
wavelengths longward of 1200 A, to iZAiZS/SOLSTICE 
solar flux measurements appro priate to the solar acti vity 
at the time of our observation (jRottman et al.ll2001l ). 

4.1.3. C/2000 WMi (LINEAR) 

The CO emission detected by STIS in the observation 
of comet C/2000 WMi (LINEAR) was too weak to al- 
low us to repeat the same analysis performed in the case 
of comets 153P/Ikeya-Zhang and C/2001 Q4 (NEAT). 
Instead, we chose to integrate the STIS spectrum over 
increasing widths centered on the nucleus, in order to 
make use of the stronger signal in the center and to in- 
crease the number of contributing pixels. We started 
with a 4" wide region which was increased progressively 
up to 16". For the rotational temperature we used the 
77 K value de rived for the cold co mponent using FUSE 
observations (| Weaver et al.l 120021 ) . while the Doppler b 
parameter value of 0. 72 km s -1 was cho sen to match the 
radio observations of iBiver et all (|2006[ ). All three STIS 
observations were averaged together to obtain detectable 
CO emission features. The best-fit column densities over 
the selected regions are plotted in Figure [9] as a function 
of the integrated slit width. Fitting a native source model 
integrated over the same rectangular regions we obtain 
a CO production rate of 3.56 ± 0.2 x 10 26 molecules s _1 . 
This model is shown by a continuous line in Figure [9l 
A CO production rate of 0.44 ± 0.03% relative to water 
is obtained adopting the favored H2O production rate 
for the FUSE ob servations (8.0xl0 28 molecules s" 1 , 
I Weaver et af]|2002T ). This makes C/2000 WMi (LIN- 
EAR) the most CO-poor comet of our sample. 

4.2. HST-GHRS Observations: C/1996 B2 
(Hyakutake) 
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Fig. 9.— CO column density profile for comet C/2000 WMi 
(LINEAR) extracted from regions of increasing widths, centered on 
the comet nucleus. The continuous line corresponds to the values 
predicted by the native source model with a production rate of 
3.56 X 10 26 molecules s _1 , when integrated over similar slit sizes. 



Since the H ST-GHRS observations do not provide spa- 
tial information, we are unable to derive a CO column 
density profile for comet C/1996 B2 (Hyakutake). The 
CO production rate is obtained by comparing the CO 
column density derived by fitting our fluorescence model 
to the GHRS spectrum with the CO column density 
predicted for the GHRS aperture by the native source 
model. As model parameters for the synthetic spec- 
trum we used a rotational temper ature of 72 K g i ven b y 
the 63 xr~ 106 K dependence from lDiSanti et al. (2003), 
which is similar to the value given bv iLis et all (1997), 
and a b parameter of 2.0 k m s -1 , within the range of out- 
flow velocities measured bv I Wouterloot et ail (|1998l ) but 
sli ghlty higher than de rived from the optical line widths 
of ICombi et al.l (|1999f ). A lower b value is inconsistent 
with the total amount of absorbed solar radiation (from 
the conservation of the number of photons), suggesting 
larger turbulent motions in the l"74xl"74 area probed 
by GHRS. 

Using the CO A — X fluorescence model with self- 
absorption we obtain a first estimate for the CO col- 
umn density averaged over the GHRS slit. However, the 
predicted line ratios are not in agreement with the data 
(see the red line in Figure [TU1 upper panel). In order to 
improve the fit and better constrain the column density 
we adjust the model to account for the geometry of the 
Sun-comet-Earth system. Starting with the previous es- 
timate on the CO column we can constrain the ratio be- 
tween the line-of-sight column density and the absorbing 
column on the comet-Sun direction. Iterating this pro- 
cedure we obtain a best fit value for the line-of-sight CO 
column density of 1.45 ±0.87 x 10 16 cm~ 2 . This model is 
shown with a red line in the lower panel of Figure[l0] The 
blue model in the same figure contains contributions from 
atomic species C II, C I, O I, and S I, as labeled, which 
account for the remaining features in the spectrum. Un- 
der the assumption of a native source model with the 
same lifetime and gas outflow velocity as employed for 
the STIS observations, we obtain a CO production rate 
of 4.97 ± 0.07 x 10 28 molecules s" 1 . This represents the 
highest CO production rate relative to water from our 
sample, 20.9 ± 0.3%, using th e water production rate of 
2.38xl0 29 molecules s" 1 from lCombi et al.l (| 19981) . 
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TABLE 3 
Production rates. 



Comet Name 


Qco* 


Qco/Qh 2 o 


Other 


Qco Measurements 




(10 28 molecules s _1 ) 


(95) 


(10 8 molecules s 1 ) 


153P/Ikeya-Zhang 


1.54 ± 0.09 


7.2 ± 0.4 a 




0.73 ± 0.16 b 


C/2001 Q4 (NEAT) 


1.76 ± 0.16 


8.8 ± 0.8 C 




1.36 ± 0.40 d 


C/2000 WMi (LINEAR) 


0.036 ± 0.002 


0.44 ± 0.03 c 




0.035 ± 0.003 f 


C/1996 B2 (Hyakutake) 


4.97 ± 0.07 


20.9 ± 0.3S 




4.84 ± 0.58 h 



* The error bars are given by the 1 a interval from the x statistics. In addition, 
we estimate that systematics amount to a 15% uncertainty in the production rates for 
the STIS observations. a Water production rate 2.15xl0 29 molecules s" 1 from HST-ST1S 
observations (see text). b Using rT^ 2 ' 1 scaling from IBiver et alj J2006). c Water production 
rate 2.0xl0 29 molecules s — 1 from HST-ST1S observations (see text). d Derived in this pa- 
per from FUSE observations (§4.1.2). The cited value includes only the cold source 
component. 6 Water production rate 8.0xl0 28 molecules s _1 used for the FUSE observa- 
tions (Weaver ct al. 2002).' Based on FUSE observations (Weaver ct al. 2002), revised for 
this paper. The cited value includes only the cold source component. 8 Water production rate 
2.38X10 29 molecules s _1 from lCombi et al.l lTl998h . h Value measured on April 1.2 using JCMT 
radio telescope IjBiver et al.l I1999h . From the 4.7x lO 28 ^" 2 ' 1 molecules s _1 dependence the 
predicted value is 6.07xl0 28 molecules s — 1 . 



Figure \W\ also clearly shows the presence of several 
bands originating on the v' = 9 and 14 levels that are 
pumpe d by solar O I A1302 and H I Lyman-q, re spec- 
tively (jWolven fc Feldmanl HifflSh . iKassall (fl976l ) first 
pointed out that scattering of Lyman-a by the (14,0) 
band is comparable to, if not larger than, the direct so- 
lar scattering by all other bands of the CO Fourth Posi- 
tive system for CO column densities > 10 17 cm~ 2 . The 
(14,4) ba nd was subsequ ently identified in the spectrum 
of Venus (|Durrancdll981h . These bands are not detected 
in any of the STIS comet spectra. The high column den- 
sity in the field-of-view of comet Hyakutake makes them 
visible, albeit at low S/N, and allows for a determina- 
tion of column density independent of optical saturation 
effects. In evaluating the fluorescence efficiency of these 
bands, the overlap between the solar lines and the indi- 
vidual lines of the CO bands is a sensitive function of 
rotational temperature and heliocentric v elocity and re- 
quire s accurate profiles of the solar lines (|Lemaire et al.1 
l200l . The (14-3), (14-4), and (9-2) bands are included 
in the model and are seen to be completely consistent 
with the CO column density derived above. 

We note that the optically thick bands connecting to 
the v" = level of the ground state, namely (1-0) at 
1509.8 A, (2-0) at 1477.6 A, (3-0) at 1447.4 A, (4-0) 
at 1419.1 A and (5-0) at 1392.6 A, all labeled in Fig- 
ure 1 101 do not show a significant variation due to geo- 
metric corrections. This can be understood from the fact 
that while the correction adjusts the line-of-sight column 
density relative to the column absorbing the solar radia- 
tion, the emission in the optically thick lines will still be 
determined by the column corresponding to one optical 
depth, which depends only on temperature and b param- 
eter. The (0-0) band at 1544.5 A does not seem to follow 
this pattern due to blending with the (3-2) band at 1542.5 
A. The same lack of variation is exhibited by the bands 
belonging to optically thin progressions pumped by so- 
lar emission lines, such as (14-3) at 1315.7 A and (9-2) 
at 1377.8 A. This is a result of the simple linear scaling 
in the optically thin limit, which makes the absorbing 
column indistinguishable from the emitting column. 
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Fig. 10.— HST-GURS spectrum of comet C/1996 B2 (Hyaku- 
take), and the model spectrum of the CO A — X system for a line- 
of-sight column density of 1.45 X 10 16 cm -2 , before (upper panel) 
and after geometric corrections (lower panel). The CO model spec- 
trum is shown in red. The blue line contains contributions from 
atomic lines of C II, C I, O I, and S I, as indicated. The A — X 
bands that do not change significantly due to the geometric effects 
are labeled. 



5. DISCUSSION 

5.1. Sources of Uncertainty and Comparison with Other 
Measurements 

The derived column densities are subject to uncertain- 
ties due to the choice of model parameters and back- 
ground subtraction. A more complete model would in- 
volve a multidimensional \ 2 minimization to constrain 
simultaneously the column density, rotational temper- 
ature and Doppler b parameter. Aside from the fact 
that this approach requires rather large computational 
resources, we expect that given the quality of the data 
the resulting \ 2 3D surfaces will have rather low con- 
trast minima and the improvement in the resulting pro- 
duction rates would be negligible. More of a concern is 
the background subtraction in the H ST-STIS data. The 
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background is variable both in the spatial and spectral 
directions from one observation to another, making it 
impossible to give a comprehensive subtraction prescrip- 
tion. The optimal background subtraction is determined 
on a case-by-case basis. While the background-related 
uncertainties can lead to ~30% variations in the values 
for the column densities, we estimate that the change in 
the resulting production rates is only about 6%. These 
values are only slightly larger than the la error bars from 
the x 2 minimization. Similarly, increasing the rotational 
temperature from 70 K to 100 K results in a ^30% de- 
crease in production rate. However, the rotational tem- 
peratures relevant for our observations were derived from 
either FUSE observations or from the T rot vs. ih de- 
pendences obtained by radio measurements, and are con- 
strained to better than ±6 K. This results only in a ^7% 
variation in production rate, as the heliocentric distance 
of the comet does not vary significantly during our ob- 
servations. The absolute values of the CO production 
rate and column density are also sensitive to the STIS 
calibration pipeline, which is based on point-source stel- 
lar standards. The column densities derived from STIS 
data could be overpredicted by at most 30% due to cali- 
bration offsets. Variations in the parameters used for the 
native source model, such as the CO lifetime and outlow 
velocity, could also change the production rate by a few 
percent. Other less quantifiable uncertainties to which 
the fluorescence model is particularly sensitive are the 
oscillator strengths and the UV solar flux, especially due 
to the variable emission features at high solar activity. 
Overall, we estimate that the systematic errors in the 
production rates amount to at most 15%. 

To assess the effects of the error sources mentioned 
above, it is useful to compare our results to other mea- 
surements of the CO production rate in these comets 
from different spectral regions. The values obtained for 
comets C/2000 WMi (LINEAR) and C/1996 B2 (Hyaku- 
take) are in excellent agreement with previous measure- 
ments. For comet C/199 6 B2 (Hy a kutak e) at simi- 
lar heliocentric distances iBiver et al.l (|1999f ) find vales 
of 4.98 ± 0.09 x 10 28 molecules s" 1 (0.952 AU) and 
4.84 ± 0.58 x 10 28 molecules s^ 1 (0.894 AU). The pro- 
duction rate relative to water is compar able to previ- 
ously measured mixing ratios of 14 to 19% (|DiSanti et al.1 
120031) and 22% (IBiver et alJfl999h. For comet C/ 2000 
WMi (LINEAR) both FUSE~ JWeaver et al J 120021 and 
radio (|Biver et al.ll2006T ) observations are consistent with 
a CO mixing ratio of ~0.4% and a CO production rate 
of ^3 to 4xl0 26 molecules s^ 1 . A good agreement is 
again found when comparing the CO production rate 
derived for comet C/2001 Q4 (NEAT) with the results 
from FUSE observations (§ 4.1.2), listed in Tabled We 
note that the measurements based on FUSE observa- 
tions used only the cold component of CO (see § 4.1.2) 
in estimating the CO production rate. The cold com- 
ponent is believed to reflect the native source of CO, 
which is directly probed by STIS. The two values agree 
marginally within the error bars, and the mismatch can 
be attributed to short time variability and pointing in- 
stability for the FU SE observations. The count rates for 
the C — X band in the FUSE data show factor of two 
variations with a periodicity of ^20 hours. The FUSE 
observation overlaps with a much shorter STIS exposure, 
but the exact correlation in comet activity between the 



two datasets is hard to assess due to the different fields 
of view of the two instruments. 

For comet 153P/Ikeya-Zhang the derived CO produc- 
tion rate is about a factor of 2 highe r than the e s timate 
from the range of values given by IBiver et al.l (|2006f) 
using the r^ 2 ' 1 scaling, see Table E]) and iDiSanti et al.1 
2002). The factor of two difference may be attributed 
to the uncertainties in the model parameters (rotational 
temperature and Doppler parameter) and in the back- 
ground subtraction. However, as discussed above, the 
combined sources of error should result in less than 
15% uncertainty. On the other hand, the r^ 2A de- 
pendence was derived from two measurements, one at 
0.51 AU and the other at 1.26 AU, not excluding the 
possibility of temporal variations affecting our obser- 
vations at 0.89 AU. Moreover, the uncertainties in the 
r^ 2 ' 1 scaling have not been quantified, making our value 
for the production rate more reliable. The CO abun- 
dance relative to water is again higher than other es- 
timates. However, using the water production rate 
given by the 19.0xl0 28 r^ 4 ' molecules s _1 dependence de- 
rived bv lBiver et al.1 (120061) from H 2 Q observat ions made 
by the Odin satellite ( Lecacheux et al.ll2003| ). the CO 
mixing ratio decreases from 7.2 ± 0.4% to 5.1 ± 0.3%. 
This would place the CO relative abundance in comet 
153P/Ik eya-Zhang closer t o the value of 3.8 ± 0.6% at 
~1 AU rtBlver et al l [111) and 4.7 ± 0.8% at 0.78 AU 
(jDiSanti et al.ll2002ft . 

5.2. Application to Broad-band Imaging 

The advantage of long-slit spectroscopy with STIS re- 
sides in giving us a better understanding of the spatial 
distribution of CO and the effects of increasing optical 
depths on the observed line ratios. We show how this 
information is applicable to broad-band imaging by in- 
tegrating the brightness of the A — X bands in the STIS 
bandpass that are not contaminated by atomic lines. 
We then integrate the fluorescence model over the same 
bands deriving an equivalent g-factor (as a ratio between 
the total brightness and the column density). Given the 
specific optical depth corrections in our model, the re- 
sulting g-factor will be in fact a brightness-column den- 
sity dependence rather than a globally constant ratio. 
The radial brightness profile of the integrated bands ob- 
served by STIS can be then converted into a column 
density profile. This procedure can be visualized in Fig- 
ure[TTJ where the data for comet 153P/Ikeya-Zhang have 
been used. The brightness profile has been rebinned in 
1" bins and background subtracted. The red histogram 
represents the column density profile obtained using the 
optically thick brightness-wavelength dependence, while 
the blue histogram represents the column density pro- 
file obtained in the optically thin approximation, using a 
constant g-factor. The black line is the native source 
model integrated in 1" x 0'.'2 bins, for a production 
rate of 1.54xl0 28 molecules s _1 , as derived in § 4.1.1. 
This figure illustrates the importance of including opti- 
cal depth effects in order to predict the correct column 
density. This method has been used to interpret the 
imaging data from the HST Solar Blind Channel of the 
Advanced Camera for Surveys (ACS/SBC), and give an 
estimate of the number of molecules released from comet 
9P/Tempel 1 as a res ult of the Deep Impact encounter 
(|Feldman et al.ll2006af h The uncertainties in this type of 



The Fourth Positive System of CO in Comet Spectra 



11 



. Native source profile 

A Model with self-absorption 
/ \ Optically thin approximation 




-10 -5 5 10 

Distance (arcsec) 



Fig. 11. — Spatial profiles of the CO column density for comet 
153P/Ikeya-Zhang, derived from the integrated brightness along 
the slit in 1" bins, selecting the unblended lines from the A — X 
system. The red line uses the total brightness-column density de- 
pendence associated with the present model, while the blue line 
is obtained in the optically thin approximation. The continu- 
ous black line is the native source model for a production rate 
of 1.54xl0 28 molecules s _1 , integrated in 1" bins. 

measurements come mainly from the additional emission 
features of atomic lines in the bandpass. 

5.3. The Native Source of CO 

The ST-STIS observations allow us to constrain the 
dominant source of CO in the region probed by the 
25" x 0"2 slit. Reconstructing the surface column den- 
sity distribution from the radial profile we estimate that 
in a 10" x 10" box the native source contribution to 
the total number of molecules is as high as 80% for 
comet C/2000 WMi (LINEAR) and from 90% to 99% 
for comets C/2001 Q4 (NEAT) and 153P/Ikeya-Zhang, 
respectively. This result suggests that the native source 
dominates in the inner ~ 3000 km from the comet cen- 
ter, while the extended component becomes important 
at larger distances. However, the native CO observed 
in the coma is linked to the amount of CO in the nu- 
cleus through the outgassing mechanism, which is not 
well known. Laboratory data suggests that the abun- 
dance of CO relative to water in the ice can be a fac- 
tor of 5 to 10 l ower t han the one observed in the coma 
(|Notesco et al.lll997t IColangeli et all 120041 ). The exact 
value depends on the temperature at which the comet is 
outgassing. 

The large variation in the CO production rate relative 
to water can be understood from the large range of he- 
liocentric distances over which the comets have formed 
and then migrated to the outer parts of the Solar Sys- 
tem through gravitational interactions with the planets. 
The observed diversity among comets is thought to be 
related to the local gas and dust composition and tem- 
perature where they formed. Under the assumption that 
CO has been trapped by water ice during comet forma- 
tion, the local temperatu re can be estimated from the 
observed production rate (jNotesco et al.lll997h . The for- 
mation temperatures in our sample are expected to range 
from 50 K or less for comets C/2001 Q4, 153P/Ikeya- 
Zhang and Hyakutake, to more than 60 K for C/2000 
WMi. If the gas trapping mechanism is more efficient 
than deposition, as sugges ted by laboratory st udies of 
CO and water interactions dCollings et al.ll2"003f) . impor- 



tant enrichments for the spe cies with a higher su blima- 
tion temperature can occur (|Notesco et al.l[T997f ). This 
can explain the lack of correlation between the abun- 
dance of CO and that of other molecules w i th a similar 
subli mation temperature (jGibb et alJ 120031 : iBiver et al.l 

Hqo2). 

6. SUMMARY 

We have developed a fluorescence model for the in- 
terpretation of CO A — X emission observed in several 
recent comets by the Hubble Space Telescope employing 
the latest values for the transition wavelengths and os- 
cillator strengths. The radiative transfer approximation 
takes into account saturation effects using Voigt profiles 
for each of the ~10 5 transitions. Self-absorption is in- 
troduced using a photon mean free path approximation. 
This process is significant for column densities above 
fewxlO 14 cm -2 , encountered close to the comet nucleus. 
It is shown that the model reproduces the optically thin 
limit for lower column densities, at larger distances from 
the comet center. When the column densities are of the 
order ~10 16 cm~ 2 or above, a good fit to the data re- 
quires a distinction between the absorbing CO column in 
the Sun-comet direction and the emitting column along 
the line-of-sight. The approximations in the radiative 
transfer model are justified for the quality of the avail- 
able data. Constraining better the optical depth effects 
demands an exact treatment, which would be warranted 
at higher spectral resolution in order to predict the rela- 
tive strengths of individual lines contained in the A — X 
bands. 

For the comets observed by STIS the fit of the column 
density profile is consistent with a predominantly native 
source, with production rates ranging from 3.56 x 10 26 to 
1.76 x 10 28 molecules s _1 . The quality and the spatial 
extent of the STIS data does not allow for a detection 
of a small extended source component. We find large 
variations in the CO abundance relative to water, from 
0.4% for C/2000 WMi (LINEAR) to 21% for C/1996 B2 
(Hyakutake). This diversity among comets is still not 
well understood, as no clear tre nd emerges at the cur- 
rent s tage of comet observations (jBockelee-Morvan et alJ 
2004). In spite of the caveats discussed in § 5, the abso- 
lute values of the CO production rate and column den- 
sity are better constrained by the introduction of optical 
depth effects, and give reasonable confidence levels for 
the CO production rate and a good agreement with pre- 
vious results. Moreover, the present model proves to be a 
valuable tool in analyzing broadband imaging data, such 
as the ACS /SBC obs ervations of comet 9P/Tempel 1 
(|Feldman et al.ll2006a|) . 

This work is based on observations with the 
NASA/ESA Hubble Space Telescope obtained at the 
Space Telescope Science Institute, which is operated by 
the Association of Universities for Research in Astron- 
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